* Wealth index calculation. Trevor Croft, 06/11/2016. Angola DHS 2015. cd 'G:\Sheas work computer\hnp2a\Angola 2015'. include file='G:\Sheas work computer\hnp2a\Angola 2015\wealth_assets.sps'. missing values QH120 (99.9). execute. set tvars=both. * Macro definitions to perform some recoding automatically. * Note: Do not use an asterisk in front of a macro name to try to comment it out - it will run any way because of the way the macros are expanded. * Delete the line instead or change the spelling of the macro in the comment. * Do not use the names of macros in comments - they will also be expanded too, and likely will fail. *{Create binary variables based on condition and give label }. define dichotomize( !positional !enclose('(',')') / !positional !tokens(1) / !positional !tokens(1) ) compute !2=0. if (!1) !2=1. variable labels !2 !3. value labels !2 0 'No' 1 'Yes'. formats !2 (f1.0). !enddefine. *(Separate shared toilets from non-shared toilets). define shared_fac( !positional !tokens(1) / !positional !tokens(1) / !positional !tokens(1) ). compute !2 = 0. do if (QH110=1). + if (!1=1) !2 = 1. + compute !1 = 0. end if. variable labels !2 !3. value labels !2 0 'No' 1 'Yes'. formats !2 (f1.0). !enddefine. *{Reset missing values to 'does not have', change 2 code to 0}. define no2zero( !positional !tokens(1) ) if (missing(!1) | !1<>1) !1=0. value labels !1 0 'No' 1 'Yes'. !enddefine. *{Construct Variables}. * Principal fonte de abastecimento de água para beber. dichotomize (QH101=11) QH101_11 "Principal fonte de abastecimento de água para beber: Torneira: Dentro de casa". dichotomize (QH101=12) QH101_12 "Principal fonte de abastecimento de água para beber: Torneira: Dentro do quintal". dichotomize (QH101=13) QH101_13 "Principal fonte de abastecimento de água para beber: Torneira: Na casa do vizinho". dichotomize (QH101=14) QH101_14 "Principal fonte de abastecimento de água para beber: Torneira: Chafariz público". dichotomize (QH101=21) QH101_21 "Principal fonte de abastecimento de água para beber: Poço/cacimba: Poço protegido". dichotomize (QH101=22) QH101_22 "Principal fonte de abastecimento de água para beber: Poço/cacimba: Poço não protegido". dichotomize (QH101=23) QH101_23 "Principal fonte de abastecimento de água para beber: Poço/cacimba: Furo com bomba". dichotomize (QH101=31) QH101_31 "Principal fonte de abastecimento de água para beber: Água de nascente: Fonte protegida". dichotomize (QH101=32) QH101_32 "Principal fonte de abastecimento de água para beber: Água de nascente: Fonte não protegida". dichotomize (QH101=41) QH101_41 "Principal fonte de abastecimento de água para beber: Água da chuva / chimpacas". dichotomize (QH101=51) QH101_51 "Principal fonte de abastecimento de água para beber: Camião cisterna". dichotomize (QH101=61) QH101_61 "Principal fonte de abastecimento de água para beber: Moto (três rodas)". dichotomize (QH101=71) QH101_71 "Principal fonte de abastecimento de água para beber: Carroça com tanque pequeno". dichotomize (QH101=81) QH101_81 "Principal fonte de abastecimento de água para beber: Lago/lagoa/riacho/canal/canal de irrigação". dichotomize (QH101=91) QH101_91 "Principal fonte de abastecimento de água para beber: Água engarrafada". dichotomize (QH101=96) QH101_96 "Principal fonte de abastecimento de água para beber: Outro". * Tipo de sanitário usam os membros do agregado. dichotomize (QH109=11) QH109_11 "Tipo de sanitário usam os membros do agregado: Dentro de casa: Sanita ligada a rede pública de esgotos". dichotomize (QH109=12) QH109_12 "Tipo de sanitário usam os membros do agregado: Dentro de casa: Sanita ligada a fossa séptica". dichotomize (QH109=13) QH109_13 "Tipo de sanitário usam os membros do agregado: Dentro de casa: Sanita ligada a fossa aberta (vala ou rio)". dichotomize (QH109=14) QH109_14 "Tipo de sanitário usam os membros do agregado: Dentro de casa: Retrete/latrina ligada a rede pública de esgotos". dichotomize (QH109=15) QH109_15 "Tipo de sanitário usam os membros do agregado: Dentro de casa: Retrete/latrina ligada a fossa séptica". dichotomize (QH109=21) QH109_21 "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Sanita ligada a rede pública de esgotos". dichotomize (QH109=22) QH109_22 "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Sanita ligada a fossa séptica". dichotomize (QH109=23) QH109_23 "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Sanita ligada a fossa aberta (vala ou rio)". dichotomize (QH109=24) QH109_24 "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Retrete/latrina ligada a rede pública de esg". dichotomize (QH109=25) QH109_25 "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Retrete/latrina ligada a fossa séptica". dichotomize (QH109=26) QH109_26 "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Retrete/latrina ligada a fossa aberta (vala". dichotomize (QH109=31) QH109_31 "Tipo de sanitário usam os membros do agregado: Fora do quintal: Sanita ligada a rede pública de esgotos". dichotomize (QH109=32) QH109_32 "Tipo de sanitário usam os membros do agregado: Fora do quintal: Sanita ligada a fossa séptica". dichotomize (QH109=33) QH109_33 "Tipo de sanitário usam os membros do agregado: Fora do quintal: Sanita ligada a fossa aberta (vala ou rio)". dichotomize (QH109=34) QH109_34 "Tipo de sanitário usam os membros do agregado: Fora do quintal: Retrete/latrina ligada a rede pública de esgotos". dichotomize (QH109=35) QH109_35 "Tipo de sanitário usam os membros do agregado: Fora do quintal: Retrete/latrina ligada a fossa séptica". dichotomize (QH109=36) QH109_36 "Tipo de sanitário usam os membros do agregado: Fora do quintal: Retrete/latrina ligada a fossa aberta (vala ou rio)". dichotomize (QH109=41) QH109_41 "Tipo de sanitário usam os membros do agregado: Balde / bacio / outro recipiente". dichotomize (QH109=61) QH109_61 "Tipo de sanitário usam os membros do agregado: Nenhum sanitário / ar livre/mato". dichotomize (QH109=96) QH109_96 "Tipo de sanitário usam os membros do agregado: Outro". shared_fac QH109_11 QH109_11_sh "Tipo de sanitário usam os membros do agregado: Dentro de casa: Sanita ligada a rede pública de esgotos - shared". shared_fac QH109_12 QH109_12_sh "Tipo de sanitário usam os membros do agregado: Dentro de casa: Sanita ligada a fossa séptica - shared". shared_fac QH109_13 QH109_13_sh "Tipo de sanitário usam os membros do agregado: Dentro de casa: Sanita ligada a fossa aberta (vala ou rio) - shared". shared_fac QH109_14 QH109_14_sh "Tipo de sanitário usam os membros do agregado: Dentro de casa: Retrete/latrina ligada a rede pública de esgotos - shared". shared_fac QH109_15 QH109_15_sh "Tipo de sanitário usam os membros do agregado: Dentro de casa: Retrete/latrina ligada a fossa séptica - shared". shared_fac QH109_21 QH109_21_sh "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Sanita ligada a rede pública de esgotos - shared". shared_fac QH109_22 QH109_22_sh "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Sanita ligada a fossa séptica - shared". shared_fac QH109_23 QH109_23_sh "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Sanita ligada a fossa aberta (vala ou rio) - shared". shared_fac QH109_24 QH109_24_sh "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Retrete/latrina ligada a rede pública de esg - shared". shared_fac QH109_25 QH109_25_sh "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Retrete/latrina ligada a fossa séptica - shared". shared_fac QH109_26 QH109_26_sh "Tipo de sanitário usam os membros do agregado: Fora de casa mas dentro do quintal: Retrete/latrina ligada a fossa aberta (vala - shared". shared_fac QH109_31 QH109_31_sh "Tipo de sanitário usam os membros do agregado: Fora do quintal: Sanita ligada a rede pública de esgotos - shared". shared_fac QH109_32 QH109_32_sh "Tipo de sanitário usam os membros do agregado: Fora do quintal: Sanita ligada a fossa séptica - shared". shared_fac QH109_33 QH109_33_sh "Tipo de sanitário usam os membros do agregado: Fora do quintal: Sanita ligada a fossa aberta (vala ou rio) - shared". shared_fac QH109_34 QH109_34_sh "Tipo de sanitário usam os membros do agregado: Fora do quintal: Retrete/latrina ligada a rede pública de esgotos - shared". shared_fac QH109_35 QH109_35_sh "Tipo de sanitário usam os membros do agregado: Fora do quintal: Retrete/latrina ligada a fossa séptica - shared". shared_fac QH109_36 QH109_36_sh "Tipo de sanitário usam os membros do agregado: Fora do quintal: Retrete/latrina ligada a fossa aberta (vala ou rio) - shared". shared_fac QH109_41 QH109_41_sh "Tipo de sanitário usam os membros do agregado: Balde / bacio / outro recipiente - shared". shared_fac QH109_96 QH109_96_sh "Tipo de sanitário usam os membros do agregado: Outro - shared". * Principal fonte de energia ou combustível que usa para cozinhar. dichotomize (QH113=1) QH113_1 "Principal fonte de energia ou combustível que usa para cozinhar: Electricidade". dichotomize (QH113=2) QH113_2 "Principal fonte de energia ou combustível que usa para cozinhar: Gás natural". dichotomize (QH113=3) QH113_3 "Principal fonte de energia ou combustível que usa para cozinhar: Petróleo / parafina / querosene". dichotomize (QH113=4) QH113_4 "Principal fonte de energia ou combustível que usa para cozinhar: Carvão". dichotomize (QH113=5) QH113_5 "Principal fonte de energia ou combustível que usa para cozinhar: Lenha/arbustos". dichotomize (QH113=6) QH113_6 "Principal fonte de energia ou combustível que usa para cozinhar: Palha/capim". dichotomize (QH113=7) QH113_7 "Principal fonte de energia ou combustível que usa para cozinhar: Cartão/papelão". dichotomize (QH113=8) QH113_8 "Principal fonte de energia ou combustível que usa para cozinhar: Fezes de animais". dichotomize (QH113=95) QH113_95 "Principal fonte de energia ou combustível que usa para cozinhar: Não cozinham em casa". dichotomize (QH113=96) QH113_96 "Principal fonte de energia ou combustível que usa para cozinhar: Outro". * Principal tipo de iluminação que utiliza. dichotomize (QH115=1) QH115_1 "Principal tipo de iluminação que utiliza: Electricidade da rede pública". dichotomize (QH115=2) QH115_2 "Principal tipo de iluminação que utiliza: Electricidade do gerador". dichotomize (QH115=3) QH115_3 "Principal tipo de iluminação que utiliza: Electricidade de luz solar". dichotomize (QH115=4) QH115_4 "Principal tipo de iluminação que utiliza: Candeeiro a pilha". dichotomize (QH115=5) QH115_5 "Principal tipo de iluminação que utiliza: Candeeiro a gás ou petróleo". dichotomize (QH115=6) QH115_6 "Principal tipo de iluminação que utiliza: Velas". dichotomize (QH115=7) QH115_7 "Principal tipo de iluminação que utiliza: Lanterna". dichotomize (QH115=8) QH115_8 "Principal tipo de iluminação que utiliza: Lenha". dichotomize (QH115=9) QH115_9 "Principal tipo de iluminação que utiliza: Não têm iluminação". dichotomize (QH115=96) QH115_96 "Principal tipo de iluminação que utiliza: Outro". * Vacas/bois. if (missing(QH118A) | QH117 <> 1) QH118A=0. missing values QH118A (99). compute QH118A_0 = (QH118A = 0). compute QH118A_1 = (QH118A >= 1 & QH118A <= 4). compute QH118A_2 = (QH118A >= 5 & QH118A <= 9). compute QH118A_3 = (QH118A >= 10 & QH118A <= 98). variable labels QH118A_0 'Vacas/bois: None' /QH118A_1 'Vacas/bois: 1-4' /QH118A_2 'Vacas/bois: 5-9' /QH118A_3 'Vacas/bois: 10+'. * Cabritos. if (missing(QH118B) | QH117 <> 1) QH118B=0. missing values QH118B (99). compute QH118B_0 = (QH118B = 0). compute QH118B_1 = (QH118B >= 1 & QH118B <= 4). compute QH118B_2 = (QH118B >= 5 & QH118B <= 9). compute QH118B_3 = (QH118B >= 10 & QH118B <= 98). variable labels QH118B_0 'Cabritos: None' /QH118B_1 'Cabritos: 1-4' /QH118B_2 'Cabritos: 5-9' /QH118B_3 'Cabritos: 10+'. * Porcos/leitões. if (missing(QH118C) | QH117 <> 1) QH118C=0. missing values QH118C (99). compute QH118C_0 = (QH118C = 0). compute QH118C_1 = (QH118C >= 1 & QH118C <= 4). compute QH118C_2 = (QH118C >= 5 & QH118C <= 9). compute QH118C_3 = (QH118C >= 10 & QH118C <= 98). variable labels QH118C_0 'Porcos/leitões: None' /QH118C_1 'Porcos/leitões: 1-4' /QH118C_2 'Porcos/leitões: 5-9' /QH118C_3 'Porcos/leitões: 10+'. * Ovelhas/carneiros. if (missing(QH118D) | QH117 <> 1) QH118D=0. missing values QH118D (99). compute QH118D_0 = (QH118D = 0). compute QH118D_1 = (QH118D >= 1 & QH118D <= 4). compute QH118D_2 = (QH118D >= 5 & QH118D <= 9). compute QH118D_3 = (QH118D >= 10 & QH118D <= 98). variable labels QH118D_0 'Ovelhas/carneiros: None' /QH118D_1 'Ovelhas/carneiros: 1-4' /QH118D_2 'Ovelhas/carneiros: 5-9' /QH118D_3 'Ovelhas/carneiros: 10+'. * Coelhos. if (missing(QH118E) | QH117 <> 1) QH118E=0. missing values QH118E (99). compute QH118E_0 = (QH118E = 0). compute QH118E_1 = (QH118E >= 1 & QH118E <= 4). compute QH118E_2 = (QH118E >= 5 & QH118E <= 9). compute QH118E_3 = (QH118E >= 10 & QH118E <= 98). variable labels QH118E_0 'Coelhos: None' /QH118E_1 'Coelhos: 1-4' /QH118E_2 'Coelhos: 5-9' /QH118E_3 'Coelhos: 10+'. * Galinhas/patos. if (missing(QH118F) | QH117 <> 1) QH118F=0. missing values QH118F (99). compute QH118F_0 = (QH118F = 0). compute QH118F_1 = (QH118F >= 1 & QH118F <= 4). compute QH118F_2 = (QH118F >= 5 & QH118F <= 9). compute QH118F_3 = (QH118F >= 10 & QH118F <= 98). variable labels QH118F_0 'Galinhas/patos: None' /QH118F_1 'Galinhas/patos: 1-4' /QH118F_2 'Galinhas/patos: 5-9' /QH118F_3 'Galinhas/patos: 10+'. * Hectares de terra. compute landarea = QH120. if (missing(QH120) | QH120 >= 99.8) landarea = 99.9. *if (missing(QH119) | QH119 <> 1) landarea=0. missing values landarea (99.9). * Electricidade. no2zero QH121A. * Rádio. no2zero QH121B. * Televisor. no2zero QH121C. * Telefone fixo. no2zero QH121D. * Computador. no2zero QH121E. * Geleira/arca. no2zero QH121F. * Internet. no2zero QH121G. * Telefone celular. no2zero QH122A. * Bicicleta. no2zero QH122B. * Motorizada. no2zero QH122C. * Carroça. no2zero QH122D. * Carro/camião. no2zero QH122E. * Barco a motor. no2zero QH122F. * Algum membro deste agregado familiar tem conta bancária. no2zero QH123. * Material principal do chão. dichotomize (QH142=11) QH142_11 "Material principal do chão: Piso natural: Terra batida/areia". dichotomize (QH142=21) QH142_21 "Material principal do chão: Piso rudimentário: Madeira". dichotomize (QH142=31) QH142_31 "Material principal do chão: Piso acabado: Tacos de madeira". dichotomize (QH142=32) QH142_32 "Material principal do chão: Piso acabado: Mosaicos de cerâmica". dichotomize (QH142=33) QH142_33 "Material principal do chão: Piso acabado: Cimento". dichotomize (QH142=34) QH142_34 "Material principal do chão: Piso acabado: Mármore/granito". dichotomize (QH142=96) QH142_96 "Material principal do chão: Outro". * Material principal do tecto. dichotomize (QH143=11) QH143_11 "Material principal do tecto: Tecto natural: Sem tecto". dichotomize (QH143=12) QH143_12 "Material principal do tecto: Tecto natural: Capim/palmeira". dichotomize (QH143=21) QH143_21 "Material principal do tecto: Tecto rudimentar: Palmeira / bambú". dichotomize (QH143=22) QH143_22 "Material principal do tecto: Tecto rudimentar: Madeira". dichotomize (QH143=23) QH143_23 "Material principal do tecto: Tecto rudimentar: Cartão". dichotomize (QH143=31) QH143_31 "Material principal do tecto: Tecto acabado: Chapas de zinco". dichotomize (QH143=32) QH143_32 "Material principal do tecto: Tecto acabado: Madeira". dichotomize (QH143=33) QH143_33 "Material principal do tecto: Tecto acabado: Lousalite/fibrocimento". dichotomize (QH143=34) QH143_34 "Material principal do tecto: Tecto acabado: Telha cerâmica". dichotomize (QH143=35) QH143_35 "Material principal do tecto: Tecto acabado: Placa de betão". dichotomize (QH143=36) QH143_36 "Material principal do tecto: Tecto acabado: Telha". dichotomize (QH143=96) QH143_96 "Material principal do tecto: Outro". * Material principal das paredes exteriores da casa. dichotomize (QH144=11) QH144_11 "Material principal das paredes exteriores da casa: Paredes naturais: Sem paredes". dichotomize (QH144=12) QH144_12 "Material principal das paredes exteriores da casa: Paredes naturais: Caniço/palmeira/paus/bambú". dichotomize (QH144=13) QH144_13 "Material principal das paredes exteriores da casa: Paredes naturais: Barro". dichotomize (QH144=21) QH144_21 "Material principal das paredes exteriores da casa: Paredes rudimentares: Pau-a-pique". dichotomize (QH144=22) QH144_22 "Material principal das paredes exteriores da casa: Paredes rudimentares: Pedra com barro". dichotomize (QH144=23) QH144_23 "Material principal das paredes exteriores da casa: Paredes rudimentares: Adobe". dichotomize (QH144=24) QH144_24 "Material principal das paredes exteriores da casa: Paredes rudimentares: Madeira". dichotomize (QH144=25) QH144_25 "Material principal das paredes exteriores da casa: Paredes rudimentares: Lata/cartão/papel/saco". dichotomize (QH144=26) QH144_26 "Material principal das paredes exteriores da casa: Paredes rudimentares: Zinco". dichotomize (QH144=31) QH144_31 "Material principal das paredes exteriores da casa: Paredes acabadas: Cimento". dichotomize (QH144=32) QH144_32 "Material principal das paredes exteriores da casa: Paredes acabadas: Pedra com cal/cimento". dichotomize (QH144=33) QH144_33 "Material principal das paredes exteriores da casa: Paredes acabadas: Tijolos". dichotomize (QH144=34) QH144_34 "Material principal das paredes exteriores da casa: Paredes acabadas: Blocos de cimento". dichotomize (QH144=35) QH144_35 "Material principal das paredes exteriores da casa: Paredes acabadas: Madeira". dichotomize (QH144=96) QH144_96 "Material principal das paredes exteriores da casa: Outro". *{Members per sleeping room}. if (hhusual=0) hhusual=hhslept. if (QH117>0) memsleep=trunc(hhusual/QH117). if (QH117=0) memsleep=hhusual. if (missing(QH117) or QH117>=99 or memsleep>=98) memsleep=99. variable labels memsleep 'Number of members per sleeping room'. value labels memsleep 0 'Less than 1 per room'. formats memsleep (f2.0). missing values memsleep (99). * Compute urban and rural variables coded (1/0) for filters later. compute urban=(QHTYPE = 1). compute rural=(QHTYPE = 2). variable labels urban 'Urban' / rural 'Rural'. value labels urban 1 'Urban' / rural 1 'Rural'. formats urban rural (f1.0). execute. * Check on indicator variable creation. frequencies variables=QHTYPE to domestic. * Toilet facility by shared/not shared. crosstabs /tables=QH109 by QH110. * Land area by units - if there are separate units - need to convert them to one unit. *crosstabs /tables=QH120N by QH120U. frequencies variables= QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea . * Turn off weights before all factor analysis. weight off. * Name the dataset window for the HH data for use later. dataset name assets. save outfile='assets.sav'. *** Test Factor Analysis. filter off. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /method=correlation. *** Common Factor analysis ***. ** Removed area-specific variables **. ** Agricultural animal variables excluded ** Land area excluded ** Any others ?. *** Common Factor Analysis. filter off. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL com) /method=correlation. *** Urban Factor Analysis. filter off. filter by urban. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL urb) /method=correlation. *** Rural Factor Analysis. filter off. filter by rural. execute. factor /variables QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea /missing meansub /analysis QH101_11 QH101_12 QH101_13 QH101_14 QH101_21 QH101_22 QH101_23 QH101_31 QH101_32 QH101_41 QH101_51 QH101_61 QH101_71 QH101_81 QH101_91 QH101_96 QH109_11 QH109_12 QH109_13 QH109_14 QH109_15 QH109_21 QH109_22 QH109_23 QH109_24 QH109_25 QH109_26 QH109_31 QH109_32 QH109_33 QH109_34 QH109_35 QH109_36 QH109_41 QH109_61 QH109_96 QH109_11_sh QH109_12_sh QH109_13_sh QH109_14_sh QH109_15_sh QH109_21_sh QH109_22_sh QH109_23_sh QH109_24_sh QH109_25_sh QH109_26_sh QH109_31_sh QH109_32_sh QH109_33_sh QH109_34_sh QH109_35_sh QH109_36_sh QH109_41_sh QH109_96_sh QH113_1 QH113_2 QH113_3 QH113_4 QH113_5 QH113_6 QH113_7 QH113_8 QH113_95 QH113_96 QH115_1 QH115_2 QH115_3 QH115_4 QH115_5 QH115_6 QH115_7 QH115_8 QH115_9 QH115_96 QH121A QH121B QH121C QH121D QH121E QH121F QH121G QH122A QH122B QH122C QH122D QH122E QH122F QH123 QH142_11 QH142_21 QH142_31 QH142_32 QH142_33 QH142_34 QH142_96 QH143_11 QH143_12 QH143_21 QH143_22 QH143_23 QH143_31 QH143_32 QH143_33 QH143_34 QH143_35 QH143_36 QH143_96 QH144_11 QH144_12 QH144_13 QH144_21 QH144_22 QH144_23 QH144_24 QH144_25 QH144_26 QH144_31 QH144_32 QH144_33 QH144_34 QH144_35 QH144_96 LAND memsleep QH118A_0 QH118A_1 QH118A_2 QH118A_3 QH118B_0 QH118B_1 QH118B_2 QH118B_3 QH118C_0 QH118C_1 QH118C_2 QH118C_3 QH118D_0 QH118D_1 QH118D_2 QH118D_3 QH118E_0 QH118E_1 QH118E_2 QH118E_3 QH118F_0 QH118F_1 QH118F_2 QH118F_3 landarea /print univariate initial extraction fscore /criteria factors(1) iterate(25) /extraction pc /rotation norotate /save reg(ALL rur) /method=correlation. * Label the created score variables. rename variables (com1 urb1 rur1=comscore urbscore rurscore). variable labels comscore 'Common wealth score' /urbscore 'Urban wealth score' /rurscore 'Rural wealth score'. * Add a variable used for linking later. filter off. string rowtype_ (A8). compute rowtype_ = 'EST'. * Calculate regressions. ** Area=urban. filter off. filter by urban. execute. * Declare a dataset to be written to in the regression. dataset declare urbancorv. * Run regression of comscore with urbscore. regression /missing listwise /statistics coeff outs R anova /criteria=pin(.05) pout(.10) /noorigin /dependent comscore /method=enter urbscore /outfile=corv(urbancorv). * Activate file of output from regression. dataset activate urbancorv. * Drop all rows of output except the coefficients. select if (rowtype_ = 'EST'). execute. * Delete unnecessary variables before merging. delete variables DEPVAR_ VARNAME_. * Rename variables containing the constant and the coefficient. rename variables CONST_=urbconst urbscore=urbcoeff. * Re-activate the main household data. dataset activate assets. * merge the coefficients. match files /file = * /table = urbancorv /by ROWTYPE_. execute. * Calculate regressions. ** Area=rural. filter off. filter by rural. execute. * Declare a dataset to be written to in the regression. dataset declare ruralcorv. * Run regression of comscore with rurscore. regression /missing listwise /statistics coeff outs R anova /criteria=pin(.05) pout(.10) /noorigin /dependent comscore /method=enter rurscore /outfile=corv(ruralcorv). * Activate file of output from regression. dataset activate ruralcorv. * Drop all rows of output except the coefficients. select if (rowtype_ = 'EST'). execute. * Delete unnecessary variables before merging. delete variables DEPVAR_ VARNAME_. * Rename variables containing the constant and the coefficient. rename variables CONST_=rurconst rurscore=rurcoeff. * Re-activate the main household data. dataset activate assets. * merge the coefficients. match files /file = * /table = ruralcorv /by ROWTYPE_. execute. dataset close urbancorv. dataset close ruralcorv. dataset activate assets. filter off. *** Calculate combined wealth score from Urban and Rural Scores. * Use coefficients from urban and rural regressions above. compute combscor=0. print formats combscor (f11.5). write formats combscor (f11.5). ** Urban. if (urban = 1) combscor=urbconst+urbcoeff*urbscore. ** Rural. if (rural = 1) combscor=rurconst+rurcoeff*rurscore. variable labels combscor 'Combined national wealth score'. execute. compute hhwt = QHWEIGHT/1000000. variable labels hhwt 'HH weights'. weight by hhwt. filter off. execute. frequencies variables=combscor comscore /format=notable /ntiles=5 /statistics=stddev mean /histogram normal /order=analysis. filter off. filter by urban. execute. frequencies variables=combscor urbscore /format=notable /ntiles=5 /statistics=stddev mean /histogram normal /order=analysis. filter off. filter by rural. execute. frequencies variables=combscor rurscore /format=notable /ntiles=5 /statistics=stddev mean /histogram normal /order=analysis. filter off. execute. *Calculate quintiles and scores for data file. compute hhmemwt=QHWEIGHT*hhusual/1000000. weight by hhmemwt. variable labels hhmemwt 'HH members weighting for index'. filter off. filter by urban. execute. rank variables=urbscore (A) /rank /ntiles (5) /print=yes /ties=mean. filter off. filter by rural. execute. rank variables=rurscore (A) /rank /ntiles (5) /print=yes /ties=mean. filter off. execute. rank variables=combscor (A) /rank /ntiles (5) /print=yes /ties=mean. variable labels ncombsco 'Combined wealth index' /nurbscor 'Urban wealth index' /nrurscor 'Rural wealth index'. value labels ncombsco nurbscor nrurscor 1 'Lowest' 2 'Second' 3 'Middle' 4 'Fourth' 5 'Highest'. *** Check on quintiles. frequencies variables=ncombsco nurbscor nrurscor. weight by hhwt. ctables /table=( QH101_11+QH101_12+QH101_13+QH101_14+QH101_21+QH101_22+QH101_23+QH101_31+QH101_32+QH101_41+QH101_51+QH101_61+QH101_71+QH101_81+QH101_91+QH101_96+ QH109_11+QH109_12+QH109_13+QH109_14+QH109_15+QH109_21+QH109_22+QH109_23+QH109_24+QH109_25+QH109_26+QH109_31+QH109_32+QH109_33+QH109_34+QH109_35+QH109_36+ QH109_41+QH109_61+QH109_96+ QH109_11_sh+QH109_12_sh+QH109_13_sh+QH109_14_sh+QH109_15_sh+QH109_21_sh+QH109_22_sh+QH109_23_sh+QH109_24_sh+QH109_25_sh+QH109_26_sh+QH109_31_sh+QH109_32_sh+ QH109_33_sh+QH109_34_sh+QH109_35_sh+QH109_36_sh+QH109_41_sh+QH109_96_sh+ QH113_1+QH113_2+QH113_3+QH113_4+QH113_5+QH113_6+QH113_7+QH113_8+QH113_95+QH113_96+ QH115_1+QH115_2+QH115_3+QH115_4+QH115_5+QH115_6+QH115_7+QH115_8+QH115_9+QH115_96+ QH121A+QH121B+QH121C+QH121D+QH121E+QH121F+QH121G+QH122A+QH122B+QH122C+QH122D+QH122E+QH122F+QH123+ QH142_11+QH142_21+QH142_31+QH142_32+QH142_33+QH142_34+QH142_96+ QH143_11+QH143_12+QH143_21+QH143_22+QH143_23+QH143_31+QH143_32+QH143_33+QH143_34+QH143_35+QH143_36+QH143_96+ QH144_11+QH144_12+QH144_13+QH144_21+QH144_22+QH144_23+QH144_24+QH144_25+QH144_26+QH144_31+QH144_32+QH144_33+QH144_34+QH144_35+QH144_96+ LAND+memsleep+ QH118A_0+QH118A_1+QH118A_2+QH118A_3+QH118B_0+QH118B_1+QH118B_2+QH118B_3+ QH118C_0+QH118C_1+QH118C_2+QH118C_3+QH118D_0+QH118D_1+QH118D_2+QH118D_3+ QH118E_0+QH118E_1+QH118E_2+QH118E_3+QH118F_0+QH118F_1+QH118F_2+QH118F_3+landarea )[S][Mean F8.3] by (ncombsco+nurbscor+nrurscor)[C] /slabels visible=no. graph /histogram(normal)=combscor /title= 'Distribution of Households by Wealth Scores'. frequencies variables=combscor /format=notable /ntiles=5 /statistics=stddev minimum maximum semean mean median mode skewness seskew kurtosis sekurt /order=analysis. weight off. filter off. write formats combscor urbscore rurscore (f11.5). * Save final dataset of assets. save outfile='assets.sav'. *** Write out scores file. write outfile='scores.dat' table /QHCLUST QHNUMBER combscor ncombsco urbscore nurbscor rurscore nrurscor. execute.